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Abstract 

A new integrable version of the degenerate supersymmetric t-J model is pro- 
posed. In this formulation instead of restricting single occupancy of electrons 
at each lattice site we may have up to two electrons at each site. As a require- 
ment of exact integrability the hopping interaction turns out to be correlated 
with the density of electrons in the neighboring sites. The exact solution of 

the model is obtained through the coordinate Bethe anzatz. 
PACS numbers: 75.10.Lp, 74.20-z, 71.28+d 
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Recently there has been considerable interest in studying low-dimensional electronic mod- 
els with strong correlation due to the possibility that the normal state of the two-dimensional 
novel superconductivity may share some interesting features of a one-dimensional interacting 
electron system. Exactly solvable models of strongly correlated electrons which are general- 
izations of the Hubbard and t- J models have been formulated and investigated [1 - 14] . It is 
noteworthy that some of these models [7,8,10] have superconducting behavior with dominant 
superconducting correlation functions. 

The t-J model is a lattice model on the restricted electronic Hilbert space, where the 
occurrence of two electrons on the same lattice site is forbidden. The Hamiltonian of the 
extended version of the t-J model with the generalized spin S — (N — l)/2 has the form 
[15,16,10] 

L N 

H = -* S ( C Ja C J+l,a + c /+l,a c J,a) 

j=l a =l 

N 

where Cj j0t and nj, a = Cj a Cj ia are the standard fermionic and density operators and Pi 
projects out from the Hilbert space the states with more than single occupancy at a given 
site. The model (1) is not integrable for all the values of the parameter J/t and arbitrary 
band-filling, but only at the supersymmetric point J — t. By construction the spin degrees 
of freedom of model (1) form a S77(iV)-invariant subspace. 

In this communication we introduce a new version of the supersymmetric t-J model. In 
this version instead of the operator Pi we introduce the operator P 2 which is a projection 
operator into the subspace where at most double occupancy is allowed at the sites. Differ- 
ently from (1) instead of N allowed states we have in this case (iV 2 + N + 2)/2 possible states 
at a given site j, namely 

|0>; c+j0>(a = l,2,...,A0; 

c+ a c+jO>(a^(3=l,2,...,N). (2) 

In a periodic lattice of length L the Hamiltonian of our model reads 
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L 
3=1 

H 3,k = J2(4,a C k,a + h.C.) 
a 

x {exp^j? Y ( n 3P + n k,p)\ ~ n 3,P n ka) 

+ Y [ u ( n j,* n 3,P + n k, a n k ,p) + t p (c+ a c+ p c k:a c k> p + /i.e.)] 

a</3 

~ J J2 c U c tp c j,P c kA Y ( n J,i + n Ki) ~ £ n jn n KS \ 

~ E ( C t<x C tp C ij C 3,l C k,pC j: a + h.C.) 
+ J E C ta C t,p C k,j C k,S C 3,5 C 3,l C k,P c k,a 

- J Y ( n j,anj,pn k)1 + n k , a n k ,pn ja ) 

+ J S n j,anj,pn kn n K5 . (3) 

Similarly as in the standard t-J model we have established the integrability of (1) only at a 
special point 

t p = J = [/ = £(e' ? -l),at?7 = m2,£ = ±l (4) 

It is interesting to point that in the case N = 2, differently from (1), the model (3) is 
integrable at an arbitrary parameter r\ [12 - 14]. 

The model (1) has been constructed in the following way. We consider all possible 
interactions that conserve the number of particles of each species separately and satisfy the 
constraint imposed by P 2 that not more than double occupancy is allowed at each site. 
Moreover we impose for the amplitude of the eigenf unctions of the Hamiltonian (1) in the 
sector with n a particles of species a(a = 1, 2, N) the following ansatz 

n 

a ± ; x n , oin) = Y A pi---p'n Qn II exp(ik Pj x Qj ) (5) 

P 3=1 

where Q = (Qi, Q n ) is the permutation of the N particles such that their coordinates 
satisfy 1 < xq 1 < xq 2 < ... < xq u < L. The sum is over all permutations P = (Pi, P n ) of 
integers 1, 2, n. In the case we have two particles at the same site xq 1 = xq 1+1 the ansatz 
(5) is replaced by 



...;x„,a n ) = 2^^Pi...PiP i+ i...Pn H exp^^xgj, (6) 

p j=i 

where the bar at the Z t/l and (/ + l) th position of the superscript indicates the pair position. 
Secondly we consider the model (3) in the sector where we have only two electrons. In this 
case the projector P 2 does not produce any restrictions and our problem is exactly equivalent 
to that of a model with only two species. The coefficients A p ^}'p n arising from the different 
permutation Q are connected with each other by the elements of the two-particle S-matrix 

N 

A..pfp 2 ... = X! ^apikpn kp 2 )Aw 5 p 2p l . (7) 
5,7=1 

A necessary condition of the compatibility of equations (7) is the fulfillment of the Yang- 
Baxter equations [17,18]. There are two integrable models with N = 2 and their S-matrices 
have such factorizable form. These models are the Hubbard model [I] and the correlated 
hopping model [6,7,12-14]. We may try to use as the fundamental building block of a general 
S-matrix of our model (3) the two-particle scattering matrix of these models. Here we use 
the second possibility and choose the S-matrix as in the correlated hopping model[12 - 14] 

s l ,{k l M) = \~\~ lPl \ (8) 

where the operator P i2 interchanges the species variables a.\ and a 2 and 

!cot \k, e = —I 
(9) 
— tan ±k, £ — +1 

To complete the proof of the Bethe anzatz (5-6) we must consider the eigenvalue equations 
in the sector where the total number of particles n — 3, 4. This gives a complicated system 
of equations for the coupling constants of the Hamiltonian (1). The solution of this system 
is presented in the Hamiltonian (3). The periodic boundary condition for the system on the 
finite interval (0,L) gives us the Bethe ansatz equation. In order to obtain these equations 
we must diagonalize the transfer matrix of a related inhomogeneous vertex model with non- 
intersecting strings [19 - 20]. The Bethe ansatz equations are written in terms of the charge 
rapidities A^ ) and additional iV — 1-spin rapidities X^(l = 1, . . . , N — 1) 
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k (0)^.- ~~ 11 \(°) _ \W 



t-t A V ~ A /3 ~ » _ TT A V~ A /3 ~ V 2 T-r A « ~ A /3 ~ Vf ^ 



_V ! = _ TT 

"i A^ - A? + i ji A^ - A^ 1; + i/2 P A^ - A^' + i/2 

Z = 1, N - 1, M = n, M N = 0, a = 1, M,, Af = A,, 

where rij = Mj_i — M,- is the number of particles of species j. The eigenenergies of the 
system are given in terms of the solutions {A^} and are given by 

n n 2 

E = -2 £ cos % = 2 £ (n - £ - ). (11) 

i=i j=i l A i J + 1 

The corresponding magnetization of the eigenstate is given by a pure Zeeman splitting of 
the SU (N) spin multiplet 

i N-l 

S z = -(N - l)n - ^ Ml (12) 
z i=i 

The ground state and excitations of the system are given by inserting in (11) the solutions 
of equations (10). The rapidities have, in general, complex values and their classification is 
analogous to that of the degenerate electron gas with an attractive (^-function potential [21]: 
(i) real charge rapidities, belonging to the set A^, and correspond to unpaired propagat- 
ing electrons; (ii) complex spin and charge rapidities, which correspond to bound states of 
electrons with different spin components; and (iii) strings of complex spin rapidities, rep- 
resenting spin states. A complex of m electrons (m < 0) is characterized [22] by one real 
£(m-i) rapidity and the remaining ones given by 

Ag) = C (m - 1) +<|, (13) 
I < m — 1, pi — —{m — I — 1), — (m — 1 — 3), (m — I — 1). 

These spin and charge strings form the classes (i) and (ii) above mentioned. In the class 
(iii) there is a set of rapidities {A^}(/ = 1,2, ...,N — 1), forming asymptotically strings of 
maximum size n, 



= A« + ifi/2 + 6®", (14) 
H = -(n - 1), -(n - 3), (n - 1). 

Here 5® M = 0(e~ aL ),a > 0, vanishes in the thermodynamic limit L — > oo and A® are 
real numbers. The above rapidities (13) and (14) are substituted into (10) and the resulting 
coupled equations for the real and A® are logarithmized. This procedure generated a set 
of integer quantum numbers for each set of rapidities: p^(C) for the real I = 0, N — 1 
and similarly for the corresponding "hole" functions. In the thermodynamic limit the Bethe 
ansatz equations reduce to sets of coupled linear integral equations for the density functions. 
After Fourier transforming these equations we obtain 

N-l 

pf H +P (0 M + E P (<?) Mexp[-(/ + g- Rg )|c|/2] 

= 2e- (m)M/2 coshcj/2, Z = 0,...,N- 1, (15) 



''cj = p^ ^(o;)exp(— m|cj|/2) 

oo 

+ E^n _1) H + *S +1) M - cosh(o;/2)^)(a;)] 



71=1 



m ,, n sinh[min(m, n)u;/2] , , . 

xexp(-max(m,n)M/2) 1 . .) \> ' J , 1 = 1,..., N-l. (16) 

The last set of equations holds for m = l,...,oo with a^,a^ h ,a^\ and being iden- 
tically zero, and pi tq = min(Z, q) — 5i >q in eq. (15). These equations are identical to those 
appearing in the degenerate electron gas with an attractive 5-function potential [21] and, 
apart from the driving terms, are also identical to those of the degenerate supersymmetrical 
t-J model [14, 22]. The expression for the energy however is different and is given by 
i w-i . (1 4-9) ] 

-E = 2s{ P - g / ^(0[ ( , + ( 2) y + (2 + Jj^lh (17) 

Next we discuss the ground state properties. First consider e — — 1. As in the ordinary t-J 
model the ground state is described by the bound complexes (13) with the maximum length 
N. In order to accommodate such complexes in the ground state we assume that the number 

6 



of electrons n is a multiple of N. In the thermodynamical limit, L,n — > oo we obtain, from 
(15-16), 

r -vo r°° N ~ l 1 2A 2A 

+ L ][g«'( 7 (A-A'))WA'MA' = ^(— ) + </(— ) 

+ / ]P(A)^A = T7P, P = 7 , (18) 



where 



«'(iA)= * 



/ y /2 + A 2 

The solution p(A) of these equations yields the ground-state energy per site 

1 r—vo roo O A OA 

l e —^-^L + L m— j+n— jWAm. (i 9 ) 

In order to solve numerically the integral equations (18) it is convenient to rewrite them in 
the following form 

1 f v o 1 

P(A) - ^ J_ vQ KM - A>(A')dA' = ^K 2 (A), 

/VQ 
p(A)dA = 2-p, (20) 
-vci 



where 



Ki{A) = K^e'^dp % = 1,2, 



~ . . e-l p l/ 2 sinhi(iV- l)p ~ , . e-l p l/ 2 sinhp 
KlW = simTIiV^ ; K2{P)= sinh§iV> " 

For the ground-state energy we now have 



1 rvo 

E = -2{p - $ 2 (0) + / K 2 (-A)p(A)dA}, (21) 

Li J — VQ 



where 



$ 2 (A) = ( e -^W + e - E ^ 1 \ p \)K 2 (p)e- ipA dp. (22) 

J — oo 

The energy density E/L (19) as a function of the density p = n/L is shown in Fig. 1 

for several values of N. For the sake of comparison we show in Fig. 2 the ground-state 



energy of the ordinary degenerated t-J model, for some values of N. These figures show 
that for small values of the density the ground-state energy of both models are the same. 
The same similarity is also expected for all densities p < 1 but for large values of N . On 
the other hand, for small values of N and large densities, where the site-exclusion effect is 
more important, these figures show quite distinct behavior for both models. 

Secondly consider the case of ferromagnetic coupling e = +1. In contrast with the 
ordinary t-J model the ground state is described by bound complexes where the rapidities 
have length 2. We have in this case an equal number of electrons with two different spin 
components. The situation is analogous to that of the correlated hopping model [14]. The 
energy density, 



The energy density E/L, as a function of p, obtained from the numerical solution of (23-24), 
is shown in Fig. 3. 

To summarise, we have presented an integrable (iV 2 + iV + 2)/2-state version of the 
degenerate supersymmetric t — J model where at most two electrons are allowed in a given 
site. We have solved the model by the coordinate Bethe ansatz method and derived the 
Bethe ansatz equations. Our results naturally raise interesting problems to be solved: the 
construction of the anisotropic version of the Hamiltonian (1) as well as its graduated version. 
In the latter case one should have an integrable model of interacting chains with spin S — 1/2 
and 1. The isotropic version for iV = 2 of this model has been solved recently by Frahm et al. 
[23] . Finally it is possible to construct new generalized integrable models by the introduction 
of the operator Pi(l > 2) which is a projection operator into the subspace where at most 
1-times occupancy is allowed at the sites. As a fundamental model we may use in this case 




(23) 



is now given in terms of the solution of the integral equations 




(24) 
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the generalization of the correlated hopping model where we have I distinct species with 
single particle as well as multi-particle hopping. For I = 3 we may use as the fundamental 
building block of the general scattering matrix the S-matrix of a model which was constructed 
recently in refs. [24,25]. For an arbitrary I, as a result of this construction, we may obtain 
the integrable model with the symmetry of quantum superalgebra U q [gl(m + 1\N — m)]k=i- 
Certainly it will be interesting to study the physical properties of these quantum chains 
particularly their phase diagram and critical exponents. 

We thank M. J. Martins and H. Frahm for useful discussions. This work was supported 
in part by Conselho Nacional de Desenvolvimento Cientifico e Tecnologico - CNPq - Brazil, 
and by The Russian Foundation of Fundamental Investigations under Grant No 99-02-17646. 
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Figure Captions 



Figure 1 - Ground-state energy per site E/L as a function of the density p = n/L for the 
exact integrable model (3) with t p — J — U — e — 1 and some values of N. 

Figure 2 - Ground-state energy per site E/L as a function of the density p = n/L for the 
degenerated t-J model (1) at J = t = 1, for some values of N. 

Figure 3 - Ground-state energy per site E/L as a function of the density p = n/L for the 
exact integrable model (3) with t p = J = U = e = — 1. 
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